Interferometric Measurement of Local Spin-Fluctuations in a Quantum Gas 
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The subtle interplay between quantum- 
statistics and interactions is at the origin of 
many intriguing quantum phenomena connected 
to superfluidity and quantum magnetism [lj. The 
controlled setting of ultracold quantum gases is 
well suited to study such quantum correlated 
systems |2j. Current efforts are directed towards 
the identification of their magnetic properties 
[3H5], as well as the creation and detection of 
exotic quantum phases [6— 8j. In this context, it 
has been proposed to map the spin-polarization 
of the atoms to the state of a single-mode light 
beam [9|. Here we introduce a quantum-limited 
interferometer realizing such an atom-light in- 
terface [lOj with high spatial resolution. We 
measure the probability distribution of the 
local spin-polarization in a trapped Fermi gas 
showing a reduction of spin-fluctuations by up to 
4.6(3) dB below shot-noise in weakly interacting 
Fermi gases and by 9.4(8) dB for strong inter- 
actions. We deduce the magnetic susceptibility 
as a function of temperature and discuss our 
measurements in terms of an entanglement 
witness. 

Quantum mechanics manifests itself in the correla- 
tions between the constituent parts of a physical system. 
These correlations quantify the probability of joint mea- 
surements, and are experimentally observable in the sta- 
tistical distribution of the outcomes of repeated measure- 
ments. Experiments studying density fluctuations have 
successfully demonstrated the potential of these tech- 
niques as a tool to study local thermodynamic properties 
of quantum gases [TTHT7] . In another context, interfero- 
metric methods have been used to study spin-fluctuations 
in atomic vapors, leading to the observation of entangle- 
ment and spin-squeezing [TU] HB HH] • More recently, sev- 
eral authors have proposed to apply similar techniques to 
map quantum fluctuations, generated by the many-body 
dynamics in a quantum gas [3-19], on to the optical field 
of a single mode probe beam. In a different approach, 
speckle noise originating from out-of-focus regions in off- 
resonant imaging was related to the spin-fluctuations of a 
Fermi gas [20]. In this letter, we use a shot-noise limited 
interferometer to directly measure the probability distri- 
bution of the local spin-fluctuations in a two-component 
quantum degenerate Fermi gas. 

Our interferometer is analogous to Young's double slit 
experiment. Two tightly focused beams, the probe and 
the local oscillator, are focused to separate points as 




FIG. 1: Interferometer Setup (a) Interferometer beams in the 
vicinity of the atomic cloud. While the probe passes through 
the cloud shown in grey, the local oscillator passes by the side 
of it. The beams (minimum l/e 2 -radius 1.2 fim) overlap in 
the far field giving rise to an interference pattern as shown. 
See methods for creation of the interferometer beams, (b) 
Optical setup to obtain two interference patterns, only one of 
which is affected by the atoms. Using a quarter-wave retarda- 
tion plate (A/4) and two polarizing beam splitters (PBS), the 
<7~-component of the polarization, which interacts with the 
atoms, is separated from the <r + -component. This yields two 
far-field interference patterns on one image as shown in the 
lower part of the figure, see methods. The lens adjusts the size 
of the patterns on the camera, (c) Level-scheme illustrating 
that only a~ light interacts with the atoms. a + light is far 
detuned from resonance. Transitions from both states to their 
respective excited states are nearly closed cycling transitions. 

shown in Figure [I] and overlap in the far field. Posi- 
tion and visibility of the resulting interference pattern 
are determined by changes in phase and amplitude of 
the probe beam, which passes through an atomic cloud, 
while the local oscillator does not. The analysis of the in- 
terference pattern thus allows the reconstruction of both 
quadratures of the probe beam, phase and amplitude, 
which carry information about the local properties of the 
atomic cloud. 

A probe beam passing through a mixture of 6 Li atoms 
in the lowest two hyperfine states |1) and |2) acquires a 
phase-shift given by 

, / niSi n 2 S 2 \ m 

^ = a °{l + 7+AS! + 1 + 7+461)' (1) 

with Co the resonant scattering cross-section, s the sat- 
uration parameter and rii and Si the line-of-sight in- 
tegrated density and the frequency detuning in atomic 
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FIG. 2: Interferometer Performance (a) Measured phase (•) 
and optical density (o) as a function of detuning. Maximum 
saturation of the probe beam s — 1.2, duration of the probe 
pulses 1.2 fj,s. Solid lines result from a fit to the phase data 
using the model described in the text yielding 712 00 = 3.2, 
effective saturation 0.6 and an effective linewidth that is 20% 
broader than the natural linewidth of 5.9 MHz, which we at- 
tribute to the probe laser. Errorbars show standard devia- 
tions, (b) Measured phase-variance (without atoms) 5(f) 2 (•) 
as a function of photon number in the probe beam determined 
from 100 measurements for each point. The empty circle (o) 
indicates the phase-variance for the intensity at which the 
spin-polarization measurements were made. Errorbars are an 
uncorrelated sum of statistical and systematic uncertainties. 
Expected phase-noise ( ) for quantum efficiency r\ = 0.6. 
The width of the line corresponds to 20% errors estimated 
from the uncertainty of our determination of rj. The gray lines 
indicate the square of the phase shift expected for a single 
atom fixed in space at a detuning of half the atomic linewidth 
for the indicated l/e 2 -waists of the probe beam |21j . 



linewidths for state \i), respectively. By choosing the 
detuning exactly in between the two resonances, Si — 
—S 2 = 6.4, the phase shift (f> is proportional to the line- 
of-sight integrated spin-polarization density m = ri\ — 
n 2 . A single measurement of the phase-shift yields the 
spin-polarization of the given experimental realization. 
Consequently, the full probability distribution of spin- 
polarization can be reconstructed from repeated mea- 
surements. 

To validate our procedure we show that (i) equation 



|T]) holds in the parameter regime of the experiment, and 
(ii) that the phase measurement is only limited by pho- 
ton shot-noise of the probe beam. Then, each addition- 
ally detected photon leads to a projection of the atomic 
state into a smaller subspace. To verify (i), we measure 
the frequency-dependent phase shift and optical density 
for a Fermi gas comprised of atoms in state |2) only, 
see methods for preparation and data-processing. Fig- 
ure [2^, shows the characteristic asymmetric profile of the 
phase whereas the optical-density is well-described by a 
Lorentzian. The solid lines result from fitting the phase 
data to equation with n\ — and agree with the mea- 
surement provided the probe duration is ~ 1 (MS and the 
saturation is less than ~ 10. The use of stronger (longer) 
pulses leads to a systematic shift of measured phases to 
larger values, due to the light-forces resulting from the 
strong focusing of the probe beam and scattering of pho- 
tons. To verify (ii), we measure the phase- variance as a 
function of the photon number as shown in Figure [2]d. 
From the power-law behavior of the phase-variance we 
deduce that photon shot-noise limits the sensitivity of 
the phase measurement. Indeed, the phase-noise is ex- 
pected to be given by 8(j) 2 = [32], where N is the 
number of photons and 77 = 0.6 the quantum efficiency, 
which is determined in an independent measurement. 

We first measure the distribution of the spin- 
polarization for weakly interacting Fermi gases at three 
different temperatures given by T\ = 8.5(5) ^K, T 2 = 
2.0(2) /iK and T3 = 0.58(5) /j>K, see methods for prepara- 
tion. Figure [3^ shows that the probability distributions 
of the spin-polarization have a narrower width the lower 
the temperature of the gas. The distributions exhibit no 
significant asymmetry and are well described by Gaussian 
functions, as expected from the large number of atoms in 
the probe beam, ~ 500 in each state. For weakly inter- 
acting Fermi gases, number fluctuations in each hyperfine 
state are independent, so that fluctuations of the spin- 
polarization are given by Sm 2 — Snl+Sn 2 ,. Consequently, 
with the onset of quantum-degeneracy, spin-fluctuations 
are reduced, since in each state only atoms at the Fermi 
edge contribute to the fluctuations, which is a manifes- 
tation of antibunching due to Pauli's principle [T3 l 114 ] . 
The measured variances of the spin-polarization, in or- 
der of decreasing temperature, are Sm^ = 35(6) /im~ 4 , 



27(3) /im 4 and 8r 



15(3) /im . Here, 
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the background corresponding to a standard deviation of 
±13(2) atoms in the probed volume has been subtracted, 
see methods. The column density in the probed region 
was rii = rt 2 = w coi/2 = 110 pw~ 2 for the gases at T 2 and 
T 3 and 20 % lower for the gas at Ti . From our measure- 
ment we determine a reduction of the spin-fluctuations as 
compared to a thermal gas prepared with the same col- 
umn density by 2.1(2) dB for the gas at T 2 and 4.6(3) dB 
for the gas at T 3 , as shown in Figure [SJd, which is in 
quantitative agreement with theory for a noninteracting 
Fermi gas, see methods. 
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FIG. 3: Spin-Fluctuations in Ultracold Fermi Gases (a) Normalized histograms for the measured spin-polarization for weakly 
interacting gases at Ti = 1.02(2) T F ( ), T 2 = 0.44(5) T F ( — ) and T 3 = 0.18(2) T F ( — ) as well as for a strongly interacting 
gas of molecules at T mo i = 0.16(7) T F ( — ), where T F is the Fermi temperature, (b) Reduction of spin-fluctuations compared 
to a thermal gas of equal column density as a function of T/T F . Data points are colored identical to (a). The solid line 
shows theory for a noninteracting Fermi gas. (c) Spin susceptibility per atom Xcoi/n co i as a function of temperature for the 
weakly interacting gas (blue points) and the strongly interacting gas (o). Also shown are the susceptibility for a classical gas 
(- -), which varies as 1/T and the theoretical expectation ( ) for the trap parameters of the gas at T3. (d) Spin-fluctuation 
density ASm 2 — kT\coi as a function of column-density n co \. The prediction for a thermal gas (- -) is shown along with the 
boundary separating the separable from the non-separable region. The blue line shows theory for a noninteracting Fermi gas 
at temperature T mo i = 0.36 /iK. Errorbars are an uncorrelated sum of statistical and systematic uncertainties. 



We now turn to the study of a gas with strong re- 
pulsive interactions, prepared close to a Feshbach res- 
onance at temperature T mo i = 0.36(10) /iK, see meth- 
ods. The resulting histogram is displayed in Figure [3^, 
and shows a distribution of the spin-polarization that 
is significantly narrower than for the weakly interacting 
gases. This reflects the fact that for interacting gases 
correlations are present between the different spin states. 
In particular, for strong repulsive interactions close to 
a Feshbach resonance, weakly bound molecules form [5] 
and the number of atoms in the two states will always be 
equal. Spin-fluctuations in the atom number difference 
are created at the cost of breaking molecules and are con- 
sequently suppressed. We measure Sm^^ — 5(2) /im~ 4 
at n co i/2 = 110 /im~ 2 , with the background subtracted 
as before. This corresponds to a reduction by 9.4(8) dB 
as compared to a noninteracting thermal gas. A reduc- 
tion by ~ 18dB is expected from ref. [9] for a molecular 
BEC at zero temperature. The observation of a lower 
value could be caused by pair-breaking or fluctuations of 
the probe frequency, see methods. 

The measured values for the spin-fluctuations can 
be used to determine the magnetic susceptibility \ via 
the fluctuation-dissipation theorem (FDT) [23 , pro- 
vided the probed system is in grand-canonical equilib- 
rium with its surroundings. Due to column-integration, 
the fluctuation-dissipation theorem here reads kT J x — 
kTxmi = ASm 2 , where Xcoi is the column- integrated 
magnetic susceptibility, k Boltzmann's constant and A 
the effective area of the probed column. For small vol- 
umes and at low temperatures corrections to the FDT are 
expected, because correlations between the probed sys- 
tem and its surroundings cannot be neglected [5U 123] . 



Using ref. [24. and including column-integration, we es- 
timate the corrections to be less than 10% even for our 
coldest samples. Figure [3J: shows the column-integrated 
magnetic susceptibility per particle x C oi/ n coi as a func- 
tion of temperature. It relates the spin imbalance to 
the energy needed to create it. For high temperatures, 
the susceptibility decreases inversely proportional to the 
temperature as expected. For low temperatures and with 
the onset of quantum degeneracy, the susceptibility satu- 
rates to a value depending on the trap details: the stiffcr 
the trap, the lower the susceptibility. The solid line shows 
theory for noninteracting fcrmions calculated for the trap 
parameters of the gas at T3. The susceptibility for the 
strongly interacting gas of molecules is lower than for the 
weakly interacting gas at T 3 , despite its lower tempera- 
ture and weaker trap. 

The link between the magnetic susceptibility and spin- 
fluctuations has been proposed as a macroscopic entan- 
glement witness in solid state systems [26]. We now apply 
this concept to our measurement of the spin-fluctuations. 
For this purpose we describe each particle as a two-level 
system, which together form an effective spin M [2"f51l27| . 
where (M z ) is proportional to the atom number differ- 
ence and SM 2 — A 2 8m 2 . For all separable states, the 
ratio of the spin-fluctuations to the total atom number 
is bounded from below. For a system in grand-canonical 
equilibrium described by a Hamiltonian which is invari- 
ant under rotations of the spin state, this bound is ex- 
pressed by the inequality 

ASm 2 2 

> o, (2) 

n co \ 6 

see methods and supplementary informations for details. 
Figure [3ji shows the spin-fluctuations A5m 2 as a func- 
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tion of the column density n co \. The gas at T 3 violates 
the above inequality, which is expected because the Pauli 
principle leads to non-separable states at low tempera- 
tures 28J. The notion of entanglement applied to in- 
distinguishable particles is subject of current investiga- 
tions (25] . 

For the case of the strongly interacting gas, we find 
= 0.11(4), violating the bound by more than 10 
standard deviations. The fluctuations are 2.4 dB lower 
than expected for a noninteracting gas at the same tem- 
perature, see bold line in Figure [3}i. This corresponds to 
2.4 dB of spin-squeezing following ref. [3"U] . 

Our analysis uses the following assumptions: (i) We 
create a gas with equal atom number in the two states, 
which leads to (M z ) = 0. (ii) We probe the system in 
thermal equilibrium. As a consequence, the transverse 
components have dephased to the amount permitted by 
Pauli's principle and we have (M x ) = (M y ) — 0. (iii) 
The time-evolution is described by a Hamiltonian that is 
invariant under rotations of the spin [31 j . It is therefore 
sufficient to measure only the z-component. 

The rigorous application of the entanglement wit- 
ness [251 [27] requires equal coupling of all atoms to the 
probe beam. Here, we take into account the inhomogene- 
ity of the probe beam by introducing the effective area 
A. 

In conclusion, we have measured the probability dis- 
tribution of the spin-flucutations in a trapped Fermi gas. 
We have discussed our measurement in terms of an en- 
tanglement witness. Our work constitutes a first step to- 
wards the detection of entanglement in many-body states 
in quantum gases. The detection of higher order correla- 
tions could be achieved by extracting higher moments of 
the probability distribution for the spin-polarization |32j . 

Methods: 

Generation of interferometer beams The interferom- 
eter beams are generated by applying two radio frequen- 
cies differing by 20 MHz along each of the axes of a two- 
axis acousto-optical deflector (AOD) very similar to pre- 
vious work [53]. This results in four beams in the -1/-1 
diffraction order of the AOD which, after passing through 
a high-resolution microscope objective, are arranged in a 
square. Two of these beams, probe beam and local os- 
cillator, have exactly the same frequency and form the 
interference pattern on the camera. The other two are 
detuned by ±20 MHz and their interference patterns av- 
erage out over the duration of the probe pulses. The in- 
tensities of the beams are controlled via the power in the 
individual radio frequencies so that the local oscillator is 
20 times as intense as the probe beam. Phase stability is 
ensured by deriving each radio frequency from the same 
source for both axes. The light beams are elliptically po- 
larized. Due to the birefringence of the atomic cloud in a 
magnetic field, also used in polarization-contrast imaging 
[34] . only the <7~-polarized component of the light inter- 
acts with the atoms, while the er + - component passes 



undisturbed, see Figure [TJd. The power ratio of o~ + - and 
<7~-component is about 10. 

Experimental sequence An equal mixture of 6 Li 
atoms in the two lowest hyperfine states, denoted by 
\l)=\mj=-l/2,mi=l) and |2)=|mj=-l/2, m/=0) is 
prepared similar to previous work [T5]. A second dipole 
trap with a wavelength of 767 nm and a l/e 2 -radius of 
10 /jm is then switched on over 500 ms and is used to lo- 
cally increase the total column-density. Finally the mag- 
netic field is ramped to 475 G where the scattering length 
is a = — 100 ao, with ao the Bohr radius. The final trap 
depths of the large dipole trap are 84 fiK for T\ , 19 fiK for 
T 2 and 10 fiK for T 3 corresponding to 1050 mW, 235 mW 
and 130 mW, respectively. The trap depths of the sec- 
ond dipole trap are 9 fiK, 4.5 /iK and 1.5 /xK, respectively. 
The central region of the cloud is probed interferometri- 
cally with a pulse of 1.2 (is duration at a maximum sat- 
uration of 9. Experiments are repeated 400 times. Im- 
ages of the whole cloud are taken after the interferomet- 
ric measurement to determine the total atom number as 
well as the temperature from the radial expansion after 
1 ms time of flight. Experiments that show a larger de- 
viation of the total atom number than 5% are discarded, 
amounting to 10 — 20% of the images. For the prepara- 
tion of the strongly interacting gas, we ramp directly to 
800 G, where a = 7000 a and the trap depth is lowered to 
17 mW followed by recompression to 30 mW, which cor- 
responds to a trap depth of 4.8 /iK for the molecules. The 
experiment is repeated 200 times. We estimate the final 
temperature to be one tenth of the relevant trap depth 
leading to T mo \ = 0.36(10) after recompression. For 
the preparation of a gas containing only atoms in state 
1 2), we hold the gas for 100 ms close to a p-wave Fesh- 
bach resonance at 159 G, which leads to the loss of nearly 
all particles in state |1). Subsequent evaporation of the 
remaining atoms leads to a non-degenerate gas of atoms 
in state |2). 

Theory We compare our measurements for the 
weakly interacting gases with theory for noninteracting 
fermions [35] . We determine the temperature and the 
chemical potential from fits to the time-of-flight images 
as in previous work [T3]. We then calculate the density 
distribution in the combined trap given the fitted tem- 
perature and atom number to determine Tp in the center 
of the combined trap. This assumes full thermalization 
in the presence of the second trap, which is supported 
by the measured spin-fluctuations. The knowledge of T, 
fi and the trap shape allows us calculate the mean and 
the variance of the atomic density along the line of sight 
using column-integration. 

Data processing Three images are obtained in each 
experiment, one with atoms in the interferometer and 
two without. The images are averaged along the direc- 
tion parallel to the fringe-pattern and a Fourier filter is 
applied to suppress the low spatial frequencies. For the 
determination of the mean phase, a simultaneous sinu- 
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soidal fit to the a~ -pattern on images with and without 
atoms is used to determine the phase. Free parameters 
are amplitude and wavelength of the fringe pattern, as 
well as a common phase for both patterns and a phase- 
difference for the picture with atoms. Residual mechani- 
cal motion, e.g. due to the microscopes, is compensated 
for by using the same simultaneous fit applied to the a + - 
pattern. For determination of the phase-fluctuations we 
find it advantageous to analyze the correlations between 
the er + - and the <r _ -pattern on each image. This allows 
us to exploit the similarity of the two patterns and use the 
cr + -pattern to noiselessly amplify the signal contained in 
the cr _ -pattern, analogous to homodyne techniques. The 
phase-shift due to the atoms causes a displacement of 
the zero-crossings of the correlation function in the im- 
age with atoms as compared to the image without atoms. 
See supplemental materials for further details on data 
processing. 

Determination of effective area The effective area A 
relates the spin-fluctuations to the mean atomic den- 
sity and corresponds to the area of a beam with uni- 
form intensity giving the same result. For the nearly 
non-degenerate gas at T\ = 1.022V the spin- fluctuations 
are well described by Poissonian statistics because the 
atoms are uncorrelated. The fluctuations are thus pro- 
portional to the average atom number in the probe vol- 
ume. This allows the determination of the effective area 
A = 0.97n co i/<5my = 4.9(8) ^m 2 , corresponding to an ef- 
fective waist of the probe beam of 1.8(3) \im. The factor 
0.97 accounts for the residual suppression of the spin- 
fluctuations at Ti = 1.022>. 

Background The contribution to the phase variance 
originating from photon shot-noise is determined from 
the two images without atoms using the above described 
procedure. This yields <5mb gr = 6.7(1.0) /im~ 4 and a 
standard deviation of t he atom number difference in 
the probed volume of \J A 2 5m^ a . i = 13(2). Frequency 

fluctuations of the probe can cause fluctuations in the 
measured phase. A variance of the probe frequency of 
2 MHz 2 would correspond to apparent spin-fluctuations 
of 5 /itm -4 at the column-density in our experiment and 
could contribute to the difference of our results compared 
to the theory in ref. [5] . 

Entanglement Witness Each atom realizes a two-level 
system represented by Pauli matrices o~ x ^ VtZ . We define 
the spin-polarization or magnetization of the probed re- 
gion as M XjV>z = Yli a x,y,z' where the sum extends over 
all atoms in the probe volume. For all separable states 
the inequality SAlf, + 8M% + 5M* > 2(N) is fulfilled, 
where (N) is the average number of atoms probed |2"5] . 
Using the invariance of the Hamiltonian under spin ro- 
tations and that we measure expectation values in the 
grand-canonical ensemble, which show the same symme- 
try as the Hamiltonian, this reduces to = ASm > § . 

J N n co i — 3 

See Supplementary Information for more details. 
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BEAM GENERATION 

The interferometer beams are generated using a two- 
axis acousto-optic deflector (2D AOD). Radio- frequency 
(RF) signals at two different frequencies are sent to both 
axes of the 2D AOD, resulting in the deflection of four 
beams in the -1/-1 order. Two RF sources (USRP 2) 
generate monochromatic signals, one at 47.5 MHz and 
the other at 67.5 MHz. The splitting and subsequent 
recombination depicted in Fig. fallows to tune the power 
ratio between the two frequencies, for each axis of the 2D 
AOD. 

After passing through a high-resolution microscope ob- 
jective, the beams deflected in the -1/-1 order are ar- 
ranged in a square and have a variable power, as depicted 
on Figure |4j The most (least) powerful of those beams 
has been deflected by the non-attenuated (attenuated) 
signals from each source. These two beams have the same 
frequency, differing by 115 MHz from that of the incom- 
ing laser beam, and they lead to the interference pattern 
observed in the experiment. Conversely, the two beams 
having the same power have frequencies differing by 40 
MHz, and contribute only to the background. 
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Voltage controlled 
attenuator 
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RF Splitter 
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FIG. 4: RF circuit for the generation of the interferometer 
beams. Two RF sources (USRP 2) generate signals at 47.5 
and 67.5 MHz, respectively. After splitting and selective at- 
tenuation, the two frequencies are combined and injected on 
each axis of the 2D AOM. Four beams are generated having 
a variable power ratio. The four spots in the image represent 
the four beams after passing through the microscope. The 
sizes of the spots represent the relative powers, and the colors 
represent different frequencies of the deflected light. 



DATA PROCESSING USING CORRELATIONS 



with the reference appears at the sides. The center part 
of this probe-reference correlation signal is selected (in 
order to avoid finite size effects), and the positions of 
the zero crossings are extracted by linear interpolation 
of the discrete signal, as depicted on Figure [5ji. The er- 
ror made by linear interpolation is typically two orders 
of magnitude lower than the photon shot-noise limit in 
our experiments. The position of a zero-crossing of the 
correlation function corresponds to the displacement of 
the er + -pattern with respect to the c^-pattern that is 
required to overlap the maxima of one pattern with the 
zero-crossings of the second pattern. The determination 
of the zero-crossings of the correlation function is thus 
a direct measurement of the relative displacement of the 
<r + -pattern with respect to the er~-pattern. 



Each run of the experiment yields three pictures : one 
taken in the presence of the atoms, and two pictures 
taken in the absence of atoms (hereafter named pictures 
1,2 and 3). The position of the crossings and thus the 
fringe displacements are measured on each of those pic- 
tures. The mean over the crossings yields the mean po- 
sition of the fringe pattern in each picture. 



In order to take full advantage of the similarity between 
the probe signal [a~ polarized) and the reference signal 
(<7 + polarized), we process the data using a method in- 
spired by heterodyne detection. Figure [l] illustrates the 
different steps of this data processing algorithm. 

The two parts of the experimental images contain the 
probe and reference signals (Figure^). A line sum along 
the direction of the fringes is first computed, yielding a 
one dimensional fringe pattern signal (Figure [5]d) . The 
right part of the signal is scaled so that the probe and ref- 
erence have roughly the same intensity. A filter in Fourier 
space is then applied to the full scaled signal, conserving 
only the Fourier components around the fringe spacing. 
The precise shape of the filter does not influence the ob- 
tained results, provided the low frequency components 
(containing the envelope of the two fringe patterns and 
the background) are removed. Figure [5p presents a typi- 
cal signal obtained after those processing steps. 

We then compute the autocorrelation function of the 
processed fringe pattern. A fixed spacing is inserted be- 
tween the probe and reference signal in the processed 
fringe pattern before the correlation function is com- 
puted. We do this so that when computing the autocor- 
relation function of the picture, the cr + - a~ correlation 
is isolated from the other contributions (er + - a + and a~ 
- er~ ) Figure [SJi presents a typical autocorrelation sig- 
nal. The autocorrelation function displays two separated 
parts. The sum of the correlation functions of the refer- 
ence with itself and the probe with itself appears at the 
center. Conversely, the correlation function of the probe 



To obtain the distribution of the fringe displacements, 
the experiments are repeated up to 400 times, over about 
2 hours. Over this period, slow drifts of the phase on 
the "atoms" pictures, of up to 5 degrees, are observed, 
most probably due to temperature variations in the en- 
vironment. In order to compensate those drifts, a sliding 
average (over 15 runs) of the positions of the crossings on 
picture 3 is computed. We take this as a measurement 
of the drift, and subtract this averaged signal from the 
positions measured on pictures 1 and 2. This subtraction 
operation amounts to measuring the displacement of each 
crossing of picture 1 to the corresponding crossing of pic- 
ture 3, corrected for long term deviations. Taking the 
average of this corrected quantity over the crossings (i.e. 
averaging the position measured for all the pattern), we 
obtain the relative displacement of the fringes on pictures 
1 and 2 compared to 3. The ratio of this mean displace- 
ment to the period of the fringes yields the phase shift 
observed on pictures 1 and 2 compared to 3. 

We now have two phase shifts measured for each run 
of the experiment, with and without atoms. Thus, we 
obtain the statistical distribution of the phase shifts in 
the presence of the atoms from all the shifts of picture 
1, together with the distribution of shifts on pictures 2, 
taken in exactly the same conditions but without atoms. 
To measure the phase variance due to the atoms, we sub- 
tract the phase variance on picture 2 (background vari- 
ance in the text) to the phase variance of the picture 
1. 
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Zero crossings 



FIG. 5: Data processing using the correlation method, a) Raw picture observed on the camera, with the reference signal on 
the left (cr + ) and the probe signal on the right (cr - ). The power ratio of the two in 20. b) The signals from the picture is 
accumulated in the direction of the fringes, yielding a one-dimensional fringe pattern. The visibility of the fringes reflects the 
ratio of probe to local oscillator, c) The probe signal on the right side has been scaled up, and the full fringe pattern is spatially 
filtered leaving only the interference contribution, d) shows the autocorrelation function of the filtered fringe, where a fixed 
spacing between the two probe and the reference has been inserted. Thus, a + - a~ correlations appears separated (at the both 
sides) from other contributions (<r + - a + and <r~ - a~ correlations). A region centered on the left side of the correlation signal, 
containing the a + - a~ correlations, is extracted. In this window, the positions of the zero-crossings of the correlation function 
are measured (see text) . 



ENTANGLEMENT WITNESS 



given direction are bounded by 



In reference [25] it is shown that the magnetic sus- 
ceptibility is a macroscopic entanglement witness. If the 
magnetic susceptibility is lower than a certain bound, the 
state describing the thermal state of the bulk system can- 
not be a convex combination of product states (separable 
state), which means that the state must contain entan- 
glement. The arguments in ref. [26] are based on the 
fluctuation-dissipation theorem which links the magnetic 
susceptibility to the fluctuations of the spins that are the 
constituents of the system. The experimental signal is 
proportional to the difference in atomic density in each 
of the two spin states. Let us define the collective spin 
describing the magnetization as 



(3) 



where <7j are Pauli matrices, 
in ref. [26] leads to 



The calculation described 



SMl + SM 2 + SM 2 > 2N. 



(4) 



Assuming that the fluctuations are the same along all 
three axes, as we will argue below, the fluctuations in a 



SM 2 > -N. 
3 



(5) 



In the experiment we have access to the column- 
integrated total density n co \ — N / A and the column- 
integrated fluctuation density 5m 2 = SM 2 /A 2 , where A 
is the effective area of the probed region. The effective 
area A is determined experimentally and takes into ac- 
count the intensity profile of the probe beam. Using this 
in ^ we get 



ASm 2 2 
> -• 

"col 3 



(6) 



The measurement of a value lower than 2/3 means that 
the state cannot be separable. 
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Invariance under Spin Rotations 



The Hamiltonian describing a two- component Fermi 
jas is given by [55] : 



h 2 v 2 

2m 



+ V ext (r)-p a *(r% 



+ J J drdr'V{r- rO*(f)j*(r 7 )j*(P)^(r) f . (7) 



Here the field operators for the two spin states denoted by 
(j G {t, i} obey the fermionic anti-commutation relation 
{^(rjo-, ^'(r')^.,} = (5 CT:(T '<5(r — r'), Vis the interaction po- 
tential, V ext the trapping potential and p a is the chemical 
potential for the two components. At the energy scale 
of our experiment, the interaction term in the Hamil- 
tonian describes the s-wave scattering of two particles, 
the strength of which is characterized by the scattering 
length a. The wavefunction describing the relative posi- 
tion of the two particles is symmetric under the exchange 
of particles. As a result the spin wavefunction must be 
antisymmetric (singlet) and is thus invariant under rota- 
tions of the spin. For a balanced gas, fif = fx^, the first 
part of the Hamiltonian is invariant under rotations of the 
spin state. The external trapping potential for the two 
components differs by the energy of the hyperfine split- 
ting. However, since no transitions between the states 
occur on experimentally relevant timescales and since we 
prepare a balanced gas, this energy difference is irrele- 
vant for the dynamics of the system. Since we prepare 
the system in a thermal state, the transverse components 
of the collective spin M have completely decohered. The 
invariance under rotations of the spin relies both on the 
structure of the Hamiltonian and on the preparation of 
a balanced, thermal state [31] 



Entanglement Witness in Grand-Canonical 
Ensemble 



Symmetry for grand- canonical expectation values In 
our experiment we measure the spin-fluctuations of a 
small subsystem of the whole atomic cloud. This sub- 
system is in thermodynamic equilibrium with the rest of 
the system. In this situation the expectation values in 
the grand-canonical ensemble of a given observable have 
the same symmetry as the Hamiltonian, which might not 
be the case for the ground state. For the specific case 
considered we have (M z ) = (UMzW), where U is an ar- 
bitrary rotation of the spin state. The state of the system 



is given by the density matrix p = -^e P H . We then have 

(UM Z U^) = TrpUMJJ^ 
1 



Z 



-Tre~ pH UM z U^ 



X -TrUe- m M z U^ {[H, U] = 0) 



Z 



\Tre~ m M z XJ^U (trace cyclic) 



Z 

-Tre^ H M z 
Z 

(M z ). 



(8) 



As an example, consider the situation of a ferromagnetic 
ground state where domains of spins pointing in a cer- 
tain direction would exist. In each realization the experi- 
ment the symmetry is broken, but the direction in which 
the spins are pointing is different. Every direction must 
be equally probable, which leads to the situation that 
the ensemble averages have the same symmetry than the 
Hamiltonian. 

Non-fixed particle number The derivation of the in- 
equality [4] assumes a fixed number of particles N and 
the authors of ref. [26 work in the canonical ensemble. 
However, the argument carries through for measurements 
in the grand-canonical ensemble where N is replaced by 
its expectation value (N). For this we express the density 
matrix in the grand-canonical ensemble as 



P 



w N p N , 

N 



(9) 



where pat is the density matrix in the subspace of fixed 
particle number N with the weight factors w^- As be- 
fore, we write (•) = Trp- for the grand-canonical expec- 
tation value and similarly we define (-)n = Trp^- as the 
expectation value in the iV-particle subspace. Similarly, 
8- and <5/y denote the variance in the grand-canonical en- 
semble and with fixed particle number TV, respectively. 
We now calculate the variance of the operator M z 



SM. 



J-2 



{{M z -{M z )f) 
^w N Trp N (M z - (M z )) 2 



N 



N 



Y^w N Trp N (M z - (M Z ) N - ((M z ) - (M z ) N )f 

(M z ) N f 

(M Z ) N ) 



^w N Trp N (M z 

N 

-2(M Z 

+ {(M Z ) - (M z/N) 

w N 8 N M 2 z + wn{{M z 



(M Z ) N )((M Z 

\2l 



(M z ) 



N 



N 



A'* 



> 



(10) 



N 
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In the last step we have used the known result from [5] 
The last term in the last line is small, because the width 
of the weights wn is centered around N — (N) with a 
width y/ (N) and because (M z ) ~ (M z )( N y More im- 
portantly, the last term is positive and if the fluctuations 
per particle are measured to be smaller than 2/3 the state 
cannot be separable. 



